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ABSTRACT 

We present our analysis of a fully sampled, high resolution dust extinction map of the Barnard 
59 complex in the Pipe Nebula. The map was constructed with the infrared color excess technique 
applied to a photometric catalog that combines data from both ground and space based observations. 
The map resolves for the first time the high density center of the main core in the complex, that is 
associated with the formation of a small cluster of stars. We found that the central core in Barnard 
59 shows an unexpected lack of significant substructure consisting of only two significant fragments. 
Overall, the material appears to be consistent with being a single, large core with a density profile 
that can be well fit by a King model. A series of NH3 pointed observations towards the high column 
density center of the core appear to show that the core is still thermally dominated, with sub-sonic 
non-thermal motions. The stars in the cluster could be providing feedback to support the core against 
collapse, but the relatively narrow radio lines suggest that an additional source of support, for example 
a magnetic field, may be required to stabilize the core. Outside the central core our observations 
reveal the structure of peripheral cores and resolve an extended filament into a handful of significant 
substructures whose spacing and masses appear to be consistent with Jeans fragmentation. 
Subject headings: infrared: ISM — ISM: dust, extinction — ISM: globules — stars: formation 



1. INTRODUCTION 

Knowledge about the initial conditions of star forma- 
tion is scarce: in a majority of the observable clouds, 
star formation has been active for relatively long peri- 
ods of time, and the clouds have evolved from their early 
stages. Examples of the few observable clouds that dis- 
play little or no star formation acti vity and may be very 
young, are the Madd alena Cloud ([M addalena fc Thad- 
deus|[l985), Lup us 5 dT achihara et aT 



199 6J, portions oT 



the Aquila Rift QKawamura et al.||1999| |Prato L || 2006 1 



and the Pipe Nebula. The latter, with a mass of 10 M , 
is particularly interestin g because it is one of the nearest 
clouds (d = 130±H pc; |Lombardi, Alves fc Lada||2006| 
that we can observe with a great level ot detail, tor these 
reasons the Pipe has recently become a prime target for 
t he study of the early stages o f star formation. 
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J (hereafter LAL06) 
Micron All Sky Sur- 



the cloud at a resolution of 60". Using this map, |Alves et 
al~1(|2007|) and|Rathborne et al.|p009| (hereafter RLA09) 



revealed and confirmed, respectively, that the Pipe Neb- 
ula hosts over 130 dense cores. According to those stud- 
ies, the mass spectrum of the Pipe Nebula cores appears 
to be very similar in shape to the standard initial mass 
function only scaled to higher mass by a factor of 3. This 
result appears to suggest that the functional form of the 
IMF may be imprinted in the primordial fragmentation 
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sure estimates for a large fraction of the cores. Analysis 
of combined data indicated that the Pipe Nebula cores 
are thermally supported, gravitationally stratified and 
pressure confined entities whi ch are possibly t he prod- 
uct of thermal fragmentation ( |Lada et al.|2008 [). A deep 
Spitzer Space telescope (SST) mid-infrared (MIPS) sur- 
vey of the Pipe cloud (Forbrich 2009, in preparation) 
show that only a handful of the more than 130 dense 
cores harbor young stellar objects, demonstrating that 
the vast majority of these cores are starless objects 

One core that stands out from the set is Barnard 59 
(B59), the largest and most massive core in the cloud. 
B59 is also one of only a few cores in the nebula with ev- 
idence for active star formation, and the only one form- 
ing a group of young stellar objects. Recent SST mid- 
infrared (IRAC) observations revealed that B59 is form- 
ing a small star cluster with 1 4 potential member s (out 
of at least 20 YSO candidates; |Brooke et alj|2007l here- 
after BHB07). Most likely the cluster is forming only 
low mass (< 2 Mq) stars. The current star forming ef- 
ficiency in the cluster would be relatively high, possibly 
above ~ 20%. 

Currently available maps of B59 are insufficient for de- 
tailed study of the core structure. The arcminute res- 
olution extinction map of LAL06 could not resolve the 
central core that hosts the cluster because the 2MASS 
observations are too shallow to penetrate the highest ex- 
tinctions in it s central region. Currently available radio 
maps of B59 ( |Onishi et aL]|1999| |Lohr et al.||2007[ ), also 
lack the resolution to resolve the internal structure of the 
core. With the goal of making a higher resolution, fully 
sampled extinction map, we combined Spitzer observa- 
tions — that penetrate the very dense center of the core 
up to a maximum column density of Ay ~100 mag — 
with near-infrared (NIR), high resolution observations 
obtained at 4 and 8m class telescopes of the European 
Southern Observatory (ESO). Our combined catalog al- 
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lowed us to determine the near to mid-infrared (MIR; 
1 to 8 /im) extinction law towards B59 to an unprece- 
dent ed depth of almost 60 visua l magnitudes of extinc- 
tion 



TABLE 1 

Color Selection for Hybrid NICE A v Evaluations 
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(I2007J). We found that the re 
sultant extinction law: 1) is approximately linear down 
to high levels of reddening, i.e., does not change signifi- 
cantly with depth; 2) departs (flattens) only moderately 
from the extinction law calculat ed with similar methods 
towards other less de nse regions (Indcbct ouw et al.|2005 
Flaherty et al.|[2~007[). These results suggest that we can 



reliably measure high dust column densities in B59 using 
a single extinction law. 

The exercise of |Roman-Zuniga et al.| gave us the confi- 
dence to construct a deep extinction map with the com- 
bined E SO-SST catalog using the standard color excess 
method ( |Lombardi fc Alves||2001 ) . The addition of mid- 
infrared colors adds the necessary power to resolve the 
central part of the region, obscured in the NIR. We il- 
lustrate this in Figure [I] where we compare measure- 
ments of dust extinction towards B59 using near-infrared 
(E(H — K)) color excess alone, and using "hybrid" 
NIR+MIR or "pure" MIR color excesses {E{m x . - K), 
E(m Xi - m Aj ); A = [3.6], [4.5], [5.8], [8.0] fim). 

The main goal of this study is to construct a deep ex- 
tinction map of Barnard 59, that resolves its central core, 
and use it as a base to further investigate the properties 
of this isolated cluster forming region. The high depth 
and resolution of the observations used in this study al- 
low us to investigate the morphology and basic physical 
characteristics of the central core and its peripheral sub- 
structures to an unprecedented level. 



Barnard 59: ESO/IRAC common region 
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Fig. 1. — Distribution of visual extinction towards Barnard 59. 
The solid histogram outlined by a dotted line represents extinc- 
tion values calculated with ESO NIR color excess, while the solid 
line histogram represents extinction values calculated with color 
excess from a combined ESO+IRAC catalog. The two vertical 
lines represent limits used in previous studies: Brooke et al. 2007 
(2MASS+IRAC) and Roman-Zuniga et al. 2007 (ESO+IRAC). 



The present paper is organized as follows: in section 
^2] we briefly describe the characteristics of the various 
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sets of data we use. In Sj3] we describe the technique 
used to construct a deep, infrared dust extinction map 
of the Barnard 59 complex. In <|4]we describe the analysis 
performed on the map to identify significant extinction 
features in it. In section Sj5]we present a discussion on 
the physical properties of significant features, their as- 
sociation with previously identified cloud cores and the 
implications for the process of fragmentation in the B59 
complex. Finally, in Sj6] we make a brief summary of our 
findings. 

2. OBSERVATIONS 

Our investigation was carried out with a combination 
of ground and space based observations. The coverage of 
the different datasets is outlined in Figure [2] 

The ground based observations are part of a high reso- 
lution NIR survey of regions across the whole Pipe Neb- 
ula. The fields corresponding to the Barnard 59 area, 
which concern this study, are divided in two sets. The 
first set are NIR (H, K s ) images obtained with the Son of 
ISAAC (SOFI) camera at the New Technology Telescope 
(NTT) facilities of the European Southern Observatory 
(ESO) atop Cerro La Silla in Atacama, Chile. These 
SOFI-NTT observations were performed at two different 
epochs: the first season ran during June 2001, with 8 
fields imaged in H and K s ; the second season was held 
in June 2002, with 9 fields imaged in H and K s . In ad- 
dition, 4 fields of the second season also have coverage in 
J. The survey aimed to cover as much as possible of the 
main structures visible in early 2MASS extinction maps. 
The resolution of the SOFI NTT observations is 0.288". 

The second set of ground-based data are NIR H and 
K s images obtained with the Infrared Spectrometer And 
Array Camera (ISAAC) at the Very Large Telescope 
(VLT) facilities of ESO atop Cerro Paranal in Antofa- 
gasta, Chile. These ISAAC-VLT observations were car- 
ried out during July 2002 with integrations of 220 and 
390 sec. in K s and H, respectively. These observations 
targeted the areas of highest extinction in the region, 
specifically the nucleus of the dense, star forming core (2 
fields) , and the westernmost core of the cloud (core 03 in 
the list of RAL09, 1 field). The high resolution (0.148") 
and sensitivity of these observations allowed us to detect 
a few more highly reddened sources that could not be 
detected with SOFI. 

Data reduction is described in detail in a following pa- 
per involving the complete ESO NIR survey of the Pipe 
Nebula (Roman-Zuniga et al 2009, in preparation). In a 
nutshell, data reducti on was carried out using a pipeline 
(| Roman-Zuniga|2007| ) that combines IRAF routines and 
ESO publicly available routines to take care of partic- 
ular aspects of the VLT and NTT data processing. A 
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Fig. 2. — Common observations of B59. The area outlined in red denotes the coverage of IRAC. The areas outlined in blue indicates the 
spawn of the NTT observations. The small areas outlines in black indicate the fields observed with the V LT. The halftone and contours in 
the background indicate extinction Ay at levels of 5, 7, 9, 12, 15, 18, 20 and 22 mag, as determined by Lombardi , Alves fc Lada] | |2006] l. 
Please see the electronic version of the Astrophysical Journal for a full color version of this figure 



second pipeline (Levinc 20061 combines the S-Extractor 
algorithm (iJertin & Arnouts 19961 with IRAF routines 
to perform Point Spread Function photometry on the 
crowded fields of the Pipe Nebula, projected against the 
Galactic bulge. 

The SST dataset corresponds to the publicly available 
pipeline catalog of B59 which is part of the Cores to Disks 
(C2D) Legacy Project and includes photometry in the 
four IRAC bands, [3.6], [4.5], [5.8] and [8.0]. Details of 
th e reduction and ca talog preparation can be consulted 
in Evans et al. ( 2005[ ) and through the SST webpages. 



n addition to the infrared data, we performed a series 
of point measurements with the Green Bank 70m tele- 
scope to determine the variation of the emission of the 
(1,1) and (2,2) rotational transitions of ammonia (NH3) 
across the extension of the central core region in B59. 
The observations were carried out during October 2007 
as a complementary project of the ammonia survey of 
dense cores across the Pipe Nebula by |Rathborne et al. 
p07t - 

2.1. Catalog Merging 

The construction of the extinction map requires a cat- 
alog that combines the NIR and MIR sets. We joined 
the photometry catalogs from different observations in 
a progressive sequence: first, the photometry lists from 
single filter observations on each instrument were com- 
bined into single-field catalogs. Second, the single-field 
catalogs were merged into single-instrument catalogs. 
Third, we constructed a complete ESO photometry list, 
by merging the ISAAC and SOFI catalogs with a source 
matching algorithm that preferentially listed ISAAC over 



SOFI detections. Positional matching between NTT 
and VLT was calculated within a tolerance radius of 
0.24", i.e. twice the average uncertainty of the ESO- 
NTT astrometric solutions. The final ESO catalog con- 
tained 59721 sources. Fourth, the ESO final list was 
merged to a 2MASS "bed" catalog obtained from the 
All Sky Data Release. The 2MASS bed catalog cov- 
ers a rectangular area from (a, 5) = (258.30,-27.70) 
to (a, 6) — (257.55,-27.55), containing 14099 sources 
after separating those fainter than the listed sensitivity 
limits and those with photometric quality flags 'U' and 
'X', which indicate upper-limits and defective observa- 
tions, respectively. In order to merge our ESO catalog 
with the 2MASS bed catalog, we determined positional 
matches in the area within a tolerance radius of 0.36" 
(three times the average uncertainty of the ESO-NTT 
astrometric solutions) and rejected the 2MASS match- 
ing sources except in a minority of cases when matching 
ESO sources were saturated. The total number of sources 
in the ESO+2MASS catalog is 71216. However, we lim- 
ited our analysis only to those sources with photometric 
uncertainty smaller than 0.1 mag in at least one band 
and no infrared excess (either from young source candi- 
dates in the cluster o r field OH IR stars, abundant in th e 
region, as shown by Lombardi, Alves & Lada (2006)). 
Our final NIR working catalog contains a total of 55257 
sources. 

Finally, we combined the NIR and the MIR data. In 
order to do this, we matched positions of sources in the 
ESO-2MASS catalog with those in the B59 STT C2D 
public catalog (5400 sources) , using a positional tolerance 
of 0.36"; this way, we found a total of 3683 NIR/MIR 
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Fig. 3. — Extinction map of the Barnard 59 region at a spatial resolution of 20". The map is presented as a contour level image with 
linear increments of 1.0 mag. The black solid line contours mark levels of extinction from Ay =6.0 to 20.0 mag in 2.0 mag increments, and 
from Ay=20.0 to 35.0 mag in 5.0 mag increments; the white contours mark levels of extinction at 40, 45, 60 and 80 mag. At this spatial 
resolution, the maximum Ay in a pixel is approximately 95 mag, however, measurements on some individual reach up to 100 mag. 



counterparts with photometric uncertainty smaller than 
0.1 mag in at least one band. However, during the 
construction of the extinction map, the ESO-2MASS- 
SST common area catalog was analyzed separately and 
merged only partially: we only used SST sources within 
an enclosed region where the density of ESO sources 
drops dramatically as one reaches the center of the cloud. 
Within this NIR obscured region, we used a combination 
of hybrid NIR+MIR and pure MIR colors for sources 
whose H or K s detections were null or had unacceptable 
photometric quality (see section [31 below) . We also re- 
jected IRAC sources which had C2D flags indicating YSO 
candidates, stars with dusty envelopes, and stars with 
flat, rising or cuspy spectral energy distributions (SED). 
In addition, from the set of objects with available [5.8] 
and [8.0] photometry, we only kept those with brightness 
within 9.0 and 14 mag, and with colors [5.8] — [8.0] < 0.5 
mag, a restriction that helps to reject candidate class III 
sources (their intrinsic values are different from those of 
background giants) and background galaxies. 

3. AN EXTINCTION MAP OF BARNARD 59 

The technique used to construct our extinction map is 
a combination of the Near Infr ared Color Excess (NICE) 
method originally presented by l Lada et al. ( 1994) a nd th e 
revised method (NICER) of |Lombardi fe Alvesj fl2001[ ). 
The NICER algorithm requires a source to be observed 
in 3 bands -in our case J, H, and K s in order to de- 
termine the extinction from two colors, otherwise it uses 
only two bands and one color, as in NICE. We do not 



have ESO J band observations available for a majority 
of the ESO fields; particularly all fields covering the cen- 
tral region. Therefore, for a majority of ESO sources, 
extinction was determined from H — K s . Also, when 
IRAC colors were added in the highest column density 
areas, we found it was more robust to use the color that 
was available for the largest number of sources possible 
at a given extinction (see below), rather than choosing 
two fixed colors from a pool with variable uniformity, 
hence the use of NICE (i.e. extinction calculated from 
one color at a time) for MIR and hybrid NIR+MIR col- 
ors. Overall, for sources using the NICE method, we 
used the following general formulation: 

A v = C el xE(\i-\ 2 ) = C e ;x[(TO Al -mA 2 ) 6 S -(m Al -mA 2 )o] 

(1) 

where m\ 1 and m\ 2 represent the observed magnitude of 
the source in any two infrared bands. The intrinsic color, 
{m\ 1 — m\ 2 )o, is calculated in a control field, in our case 
a low column density area located approximately 2° west 
of B59 at the same galactic latitude. We also have ESO 
observations and a C2D catalog for this field. C e i is a 
constant derived from the adopted extinction law, in the 
form: 



C P i — 



(2) 



A Ks \A Ks A K 

We used Ay / A k . = ■ 1 1 2 from the extinction law of |Rieke| 
fc Lebofsky| ( 1985 ), and the coefficients A\/Ak s are those 
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Fig. 4. — a A v vs Ay distribution for the extinction map of 
figure [3] Individual measurements per pixel are shown with gray 
dot symbols. The black dot symbols with error bars joint by a 
dashed line indicate the median values in 5 magnitude bins. Notice 
the spiking of the noise amplitude for Ay > 60 mag. 
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For fields where J, H and K s are available (the 2MASS 
"bed" catalog described above) Ay can be estimated 
from two colors, J — H and H — K s , by applying the 
NICER formulation of |Lombardi fc Alves| supposing 
that the extinction can estimated as Ay = a + 61 ( J — 
H) + b2(H - K) then 



A v = b x [E{J - H)) + b 2 [E(H - K)] 



(3) 



and the coefficients b\ and 62 are calculated through a 
minimum variance method. In NICER, Ay is required 
to be an unbiased estimator, i.e. b\/C e i 1 + b2/C e i 2 = 1 
and a + bl(J — H) a + b2{H - K) = 0. ' 

Ay is evaluated from the J, H and K s observations 
in most of the survey areas, where column densities were 
below ~ 4 x 10 22 cm -2 {Ay ~ 40 mag). At the dense 
central region of B59, column density rises rapidly to 
values twice as high. We delimited the highest den- 
sity region as a rectangular area running from (a, 5) = 
(257.810,-27.470) to (a, S) = (257.880,-27.390), i.e. 
slightly larger than the area observed with the VLT. In- 
side this region we used ESO+IRAC counterparts fol- 
lowing a protocol: a) if a source had reliable H and K s 
photometry, we evaluated Ay from E(H — K s ); b) if a 
source was null or unreliable in H then we evaluated Ay 
from ^([S.6] — K s ); c) if a source was null or unreliable 
in K s we evaluated Ay from i?([3.6] — [4.5]); d) finally, 
if a source was unreliable or null in [3.6] we evaluated 
Ay from _E([4.5] — [5.8]). This protocol is summarized in 
Table [T] where we listed the corresponding values of C e z 
and the intrinsic values (m\ 1 — toa 2 )o obtained from our 
control field. 

3.1. Construction of Map 

The complete extinction map, presented in Figure [3] 
was constructed using a Gaussian spatial filter with 
FWHM=20", the beam width was increased to 24" (20% 
oversampling) in the Ay > 40.0 mag region. The map 
is Nyquist sampled, yielding pixels separated by 10". 
The value of extinction at each position is the Gaussian 
weighted median of values for individual sources within a 
radius equal to 3-sigma of the Gaussian beam. The error 
in the mean extinction, a a v , was calculated as the stan- 
dard deviation of individual background column density 
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Fig. 5. — An illustrative sequence of the wavelet filter process 
on the B59 extinction map. From top to bottom: the original 
extinction map, a fifth scale wavelet filtered map, and a residual 
map resulting from the subtraction of the second from the first. 



measurements, c*, divided by the square root of the num- 
ber of sources (a*/y/N]*) within the 3-sigma circle. No 
additional smoothing was applied to the map. In the pa- 
per, a simple interpolation routine (IDL TRI-SURF) was 
used to smooth out contour jagging. The resultant map 



improves the resolu tion of the Ak map of Lombardi 
Alves & Ladal (120061) by a factor of three, allowing us 



to resolve for the first time: a) the detailed structure of 
the central core, b) the detailed structure of the thick fil- 
ament located east of the central core (hereafter eastern 
filament), c) the structures of several, peripheral cores 
located to the southeast and west of the cloud, and d) 
the faint filaments that connect the main core with the 
peripheral structures. 

The spatial resolutions of 20 and 24" for this whole 
area map were chosen so that we could assure 1) that 
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the average extinction per pixel could be calculated with 
at least two stars per pixel all the way to the peak of 
the central core; 2) that the map reveals most of the 
features visible at a threshold resolution (18"), at which 
a number of empty pixels start to appear. Empty pixels 
would especially affect the central core region, where the 
number of NIR sources starts to decrease significantly for 
Ay > 45. 

3.1.1. Map noise 

Assuming that the quality of the photometry is uni- 
form across the map, it is expected that the error in the 
measurement of extinction at a given line of sight — i.e. 
the noise per pixel — , will be inversely proportional to 
the number of sources involved. In Figure [4] we show 
a plot of the error in the mean extinction, o~a V i as a 
function of Ay. In general, at low (Ay < 10 mag) and 
moderate (10< Ay <30) extinction levels, the error is 
low, 0.11 ± 0.04 mag, and also uniform, indicating that 
any structures within those levels are well resolved at the 
chosen resolution. Above Ay ~30 mag, a a v increases 
monotonically until it reaches a maximum median value 
of about 1.5 mag per pixel within the most obscured 
regions in the central core (Ay >75 mag). At high ex- 
tinction levels, the significance of features in the map is 
determined by this monotonic increase of the noise am- 
plitude. For example, a 5a a v interval in low to moderate 
extinction regions (10 < Ay < 30 mag) increases slowly 
from 0.55 to 1.5 mag, but has to be increased from 3.0 to 
5.0 mag at the highest extinction levels (Ay > 50 mag). 
This indicates that any features in the map, especially 
peaks and valleys indicating the presence of dense gas 
cores, are only to be trusted if they can be measured 
as variations in extinction larger than the intervals just 
described. 

4. IDENTIFICATION OF DENSE CORES 

A dust extinction map is a two-dimensional (2D) struc- 
ture. Potential dense cores can be identified as local max- 
ima, forming extinction peaks that merge with the cloud 
at a local background level and have roughly elliptical 
shapes. The identification of individual cloud cores in a 
2D map would be relatively straightforward if the extinc- 
tion peaks were always well separated from each other. 
In practice, neighboring peaks may be seen in projec- 
tion close to each other and in some cases overlap, mak- 
ing it difficult to determine their boundaries at the local 
background level. In order to reduce this effect, we re- 
duced the contribution of the local background by apply- 
ing a wavelet transform filter to the map using a r outine 



developed by Benoit Vandame (see Alves et al. I . The 
routine is based on a Multi-scale Vision Model ( MVM) 



designed to reconstruct astronomical images (Rue & Bi 



jaoui 



T997J). The MVM filtering is performed by isolat- 
ing features rising above a 4. Oct threshold in the wavelet 
transform space at five spatial scales: 0.4, 0.8, 1.6, 3.2 
and 6.4'. At each scale, the pixels with values higher 
than the threshold in the wavelet transform space are 
identified and labeled, and the program constructs a se- 
ries of trees of interscale connectivity, which are then 
used to identify significant objects and to reconstruct 
the image. The resultant MVM filtered image contains 
zero or minimum background, making it easier to define 
the boundaries of individual clumps. A simple check was 



done by subtracting the filtered image from the original 
map: the resultant residual image contains only back- 
ground extinction and low level extinction from a few 
filamentary structures, mostly located east of the central 
core (see Figure [5]). 



4.1. Significant Features 

Identification of individual clumps in 2D maps (dust 
extinction and dust emission maps are common exam- 
ples) is usually done by means of a clump finding algo- 
rithm. 2D maps lack information on the velocity dis- 
tribution of the gas at each line of sight, making it dif- 
ficult to separate clumps if they overlap. Fortunately, 
the case of the Pipe is particularly benign because the 
cloud is projected across its longest axis on the plane 
of the sky. Most clump finding algorithms depend on 
a set of parameters to success at detecting peaks and 
to define the boundaries of individual features, and thus 
their reliability is ultimately subject to the correct choice 
of tho se parameters. We u sed the CLUMPFIND-2D algo- 
rithm ( Williams et al.|1994 1 , because it has only one free 
parameter, the set of contour levels used to identify ad- 
jacent regions associated with a local peak or maxima. 
The intervals used to define the contour levels are defined 
by the user and have to be chosen carefully in order to 
select only significa nt structure and avoid id entification 
of spurious featu res ( Kainulainen et al.|2009 1 . Moreover, 
as discussed in §4.3[ additional criteria are required to 
determine the significance of individual detections. 

To define the best set of contour levels, we consid- 
ered the increase of the noise amplitude as a function 
of extinction in our map. We chose intervals starting at 
Ay =0.55 mag contour level, equal to 5 times the average 
noise per pixel in low extinction (Ay < 5 mag) regions. 
For low and moderate extinction levels (1< Ay < 10 
and 10 < Ay < 25) we used fixed 5<ja v intervals set 
to 0.55 and 1.25 mag respectively. For thicker regions 
we calculated the 5cta v interval corresponding to the in- 
crease in the noise as shown in Figure [4] until reaching 
Ay=73 mag. In a first run, we kept a list of all cores 
detected with the chosen set, then we decreased the level 
values by 1 and 2 sigmas (test level dropping) in order 
to discard spurious detections. A few features were no 
longer identified after the level dropping and were con- 
sidered as spurious detection. Also, we increased the 
minimum number of pixels for region identification from 
20 to 50, to avoid confounding single pixel fluctuations 
with extinction peaks. Our final CLUMPFIND-2D run yield 
20 significant extinction peaks in the map, as shown in 
Figure 6| A list of these peak identifications is presented 
in TablelH 



4.2. Feature Properties 

For each of the extinction peaks identified, 
CLUMPFIND-2D defines an optimized boundary at 
the bottom contour level (0.55 mag, in our case). The 
size of a core or peak gets defined as the equivalent 
radius of the boundary. Mass was estimated by sum- 
ming the background corrected total extinction in 
pixels within the boundary. The conversion to mass, 
assuming a standard value for the gas to dust ratio 
Nn/Ay = 2.0 x 10 21 cm~ 2 , is given as: 
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Fig. 6. — Identification of main cloud features in Barnard 59. The map shows the location of cores identified with CLUMPFIND-2D 
on a wavelet filtered map. The 5<r intervals used to identify features in the map as shown as iso-contour levels. Please see the electronic 
version of the Astrophysical Journal for a full color version of this figure 



M core = 1.28- 10 (pj 2 {^^ 2 j^f l{ A v )M Q (4) 

Where J^fLj. /iO^v) adds the contribution of N pixels 
within the feature boundary in the wavelet filtered map, 
9 is the beam size, and D c i ou d is the distance to the cloud. 
In our case, D =130 pc and 9 = 20". 

Using the values for mass and equivalent radii, we 
made estimations of the average density, n of the features 
and, assuming a gas temperature of 10 K we calculated 
the Jeans length and the Jeans mass for each feature as: 

/ T \ 1/2 / n \- x / 2 

and 

/ T \ 3/2 / n \~ 1/2 

M ' = 1M {m) (lows) [M ° ] - (6) 

4.3. Significance of features 

The next task was to determine which of the features 
we identify correspond to cloud cores previously identi- 
fied in the study of RLA09. The core list of RLA09 for 
the same region contains only 7 cores, resulting of merg- 
ing some of the individual features they identified in the 
map of LAL06. The merging of features in RLA09 was 



done by following a well defined prescription: a feature is 
significant if its peak rises by at least 3er above the noise, 
and two features represent independent cloud cores only 
if their separation is larger than the mean Jeans length 
of the cloud (approximately 0.26 pc) or if their central 
velocities (measured from C 18 observations) differ by 
more than the 1-D isothermal sound speed at T = 10 K. 

Each of the 20 extinction peaks identified in our map, 
coincide with at least one of the 7 cores defined by 
RLA09. Therefore, many of the features in our map 
may not be independent cores but they may represent 
significant substructure within larger, independent cores. 

Our high resolution extinction map allows us to mea- 
sure improved masses and sizes of features identified with 
CLUMPFIND-2D and therefore provides improved measure- 
ments of the local Jeans length (LJL) in the region, that 
we can then compare to the peak to peak projected sep- 
arations of features which are also slightly better defined 
in our map. In what follows, we will discuss the sig- 
nificance of the features found in our map and, we will 
make a distinction between substructure peaks and inde- 
pendent cores based on the rules of RLA09. We cannot 
apply the rule of the difference of central velocities be- 
cause we do not have pointed radio observations for each 
of the peaks we identified, but we can use the projected 
peak to peak separations in our map to determine if any 
individual extinction peaks are separated by more than 
the LJL. 
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Barnard 59 NICE 2MASS + NTT+VLT+SST Dust Extinction Map 




257.90 257.85 257. E 

Right Ascension 

Fig. 7. — A close look map of the central core in B59, at the same 
resolution as the map in Figure [3] The diamond shaped symbols 
indicate the locations of YSO candidates identified by Brooke et 
al. (2007). The circle shaped symbols indicate sources identified 
as Class 0/1. Please see the electronic version of the Astrophysical 
Journal for a full color version of this figure 



4.3.1. Central Core: structure and relation to young cluster 

members 

In Figure [7] we show a close-up map of the central core 
region. We also marked the locations of the YSO candi- 
dates identified by BHB07 4 . We confirm that 13 of the 
20 YSO candidates are located within the central core, 
specifically above extinction levels of Ay = 15 mag, and 
out of these, 12 are located below Ay = 50 mag. Only 
one star, a class 0/1 source (source No. 10 in BHB07 list) 
is projected towards an extinction higher than Ay = 70 
mag. 

The central core is approximately trapezoidal in shape; 
the Northwest corner of the trapezoid appears to stand 
out from the main body at S = —27.40, defining the 
core 09c. This core is significant, with at least three 
closed contours at the 5<ta v level between Ay=25 and 
30 mag. Moreover, the projected distance from 09c to 
the midpoint between 09a and 09b, is about 0.12 pc, 
slightly longer than the LJL at the central region. Core 
09c is located next to what seems to be a cavity opening 
around (a, S) = (257.86,-27.40) -coincident with source 
13 in BHB07 list. The cavity in the cloud is likely to be 
related to the gas o utflow, evidenced by the ample lin e 
wings of 12 CO lines ( lOnishi et al.|l999||Lohr et al.|2007 1 . 
The most likely candidate responsible fo r exca vating the 
cloud is source 2M171123 ( |Riaz et al.1|2009[ ), listed as 
Source No. 11 in BHB07 list, also asso ciated with the 
millimeter emission source B59-MMS1 (Reipurth et al. 
19961. In the IRAC images, a significant scatter light 




257.85 
Right Ascension 



Fig. 8. — A close look of the central core region in the identifi- 
cation map of figure [6] The 5a intervals used to identify features 
in the map as shown as iso-contour levels. From top to bottom, 
the five top levels are Ay =73,66,60,55 and 50 mag. Please see 
the electronic version of the Astrophysical Journal for a full color 
version of this figure 



feature appears to confirm the action of the source on 
the surrounding material. 

Our map appears to show a second indentation in the 
core at (a, S) = (257.842, -27.478), just south of source 
No. 9 in BHB07 list (2M171121), which also shows a 
significant scatter light feature in NIR images a nd it 
is possibly associated with an outflow structure (Riaz 
et al.|[2009[ ). If so, the indentation in the map could 



4 Additional information as well as list of ESO deep level pho- 
tometry and a color-magnitude diagram for the cluster candidates 
is presented in a subsequent spectroscopic study of the Barnard 59 
cluster (Covey et al., in preparation) 



be an indication of another cavity like structure being 
formed by a protostellar outflow. Outflows from young 
stars have been long known to play an imp ortant role 
in the evolution of star formin g clouds (e.g. |Bally et al. 
|1999|[Matzner fc McKee|2000[ ). Fragmentation enhanced 
by outfl ows has also been pred icted by recent numerical 
studies ( |Nakamura fc Li||2007 l. 

Regarding the central region, our map reveals struc- 
tural details that could not be previously resolved in 
lower resolution maps. In Figure [8] we show a "close 
up" of the identification map with the 5a^ v intervals 
used for CLUMPFIND-2D. Starting from the top, at the 
Ay —73 mag contour, the core appears to split in two 
with source No. 10 located in the middle. The central 
split is significant at a 5cr J 4 v level at the maximum local 
noise amplitude (2.0 < {(Ja v ) < 3.0 if 65 < Ay < 80 
mag). At the next contour level, Ay =66, the contours 
close into a single feature, and no more substructure is 
found until we go down to the 40 < Ay < 50 mag range. 
There, the map shows some small closed contours, but we 
found that they do not represent significant features, as 
they cannot longer be detected if the levels are dropped 
by lcr. 

In summary, CLUMPFIND-2D suggests that the B59 core 
might be split into three fragments: two defined by the 
split in the filtered map at Ay=73 mag (listed as 09a and 
09b in the identification list of Table [2]) , and one near 
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the northern cavity (listed as 09c in the identification 
list). However, in section 5.1 we argue that the split 



near the center might not be as physically significant 
as CLUMPFIND-2D suggests, and that the central core is 
likely a single stratified core with only one independent 
new fragment (09c) identified in this high resolution map. 

4.3.2. Peripheral Cores 

Four cores are located in the periphery of the com- 
plex, corresponding to cores 3, 4, 5 and 13 in the list 
of RLA09. These objects have a morphology reminis- 
cent of Bok globules, but in our map, they appear to 
be still connected to the rest of the cloud region by thin 
(3 < Ay < 6 mag) filaments. Core 03 is clearly centrally 
condensed, with an elongated N-S distribution near the 
center; it also has a very extended "tongue" at its north- 
ern edge which makes it 40" wider at the top; its peak 
extinction is 19.9 mag in the filtered map. Core 05 has 
a peak extinction of 10.9 mag in the filtered map and 
has a toroidal li ke structure, simi lar to Globule 2 in the 
Coalsack cloud |Lada et al.||2004[ ), which CLUMPFIND-2D 
picks as a double core. We merged the two fragments, 
05a and 05b because their projected separation is shorter 
than their LJL. In the case of peaks 04a and 04b, their 
morphology in the map and optical images might suggest 
they are independent, but as their projected separation 
(0.10 pc) is about half their LJL (~0.2 pc), we cannot 
tell if they are independent. Thus, we merged features 
04a and 04b into a single core, 04ab. 

There are two highly elongated clumps at the North- 
western edge of the B59 region, both significantly 
larger in the E-W direction than in the N-S direction. 
CLUMPFIND-2D suggests these cores are each divided in 
two fragments (cores 09d,e and 07a, b). The peak-to- 
peak projected distance from 09d to 07a is about 270" 
or 0.17 pc, similar to the LJL of both cores, thus we 
confirm that peaks 09 and 07 are independent . 

In the case of peaks 07a and 07b we find that the pro- 
jected distance (0.06 pc) is shorter than the LJL, thus we 
merged into a single core, 07. Same criterion is applied 
to peaks 09d and 09e, which we merged into core 09de. 

In the case of peaks 09d and e, the analysis of RLA09 
suggested they share the same central velocity value as 
the central core, 09a-c, and for that reason they were 
merged to it. However, in our map, the projected dis- 
tance from the midpoint between 09d and 09e to the 
midpoint between 09a and 09b is 0.16 pc, which is com- 
parable to the LJL and larger than the LJL of 09a and 
09b. Thus, we consider peak 09de to be independent 
from the central core. 

In RLA09, peak llg was merged to core 11 (the east- 
ern filament) as it was found to have the same central 
velocity. However, we find that the gap between peaks 
llg and 11a in the map is very significant (in optical 
and NIR images the gap is transparent to background 
stars), and that the projected separation from peak to 
peak, 0.18 pc, is about the same as the LJL in both 
cores (0.17-0.18 pc). Thus, we suggest that llg is an 
independent core. 

The eastern filament in B59 extends east to west 
from (a, S) = (257.9377, -27.4193) to (a, S) = 
(258.1793, -27.3113). The length from tip to tip, as 
projected on the map, is approx. 0.57 pc. Its eastern- 
most tip, listed as 11a (core 11 in the list of RLA09) 



has a peak extinction of 20 mag. We identified 6 extinc- 
tion sub-peaks across the filament, which we find to be 
suggestive of significant sub-structure, as we discuss in 



(5.4 



5. A BRIEF DISCUSSION ON FRAGMENTATION IN B59 

The high spatial resolution achieved in our map should 
allow us to determine the presence of substructure re- 
lated to fragmentation in the central core of the Barnard 
59 complex. However, the map does not show compelling 
evidence of multiple fragmentation, which might be ex- 
pected from a core that has already formed more than 
a dozen stars and has been active for at least the age of 
the stars, estimated to be 2±1 Myr (Covey et al 2009, in 
preparation) . 

5.1. Fragmentation of the Central Core: A Lack of 
Substructure 



As described in (4.3.1 CLUMPFIND-2D divides the re- 



gion below the cavity in the central core in two pieces 
given that there are two closed contours at the Ay=72 
mag, the highest contour level we used. Put together, the 
fragments have a total mass of 19 Mq, which is about 10 
times as large as the critical B onnor-Ebert mass of the 
cloud, 2 M© (Lada et al. 20081 and almost 20 times as 
large as the local Jeans mass, estimated to be close to 
1 Mq. However it is not clear whether the split might 
represent the existence of independent cores: the split 
at the Ay=72 mag level is significant above the local 
noise amplitude, but as we move away from this central 
region towards lower column density values, we cannot 
find any other evidence of significant substructure. The 
total extinction above the second highest contour level 
at Ay—6Q, is 2422 mag, distributed over 33 pixels, and 
yielding a mass of approximately 2.1 Mq -that is, if we 
include the column density in those pixels all the way to 
Ay=l mag. However, if we count the extinction using a 
baseline level of A v =66 mag, the two peaks at ^4y=73 
would represent two relatively small fluctuations above 
a smooth overall density gradient in a single stratified 
massive core. In such case, the extinction above Ay—7S 
mag would be equivalent to 0.55 M Q , or about 3% of the 
total mass in the core 09ab. 

The cores 09a and 09b have a peak to peak projected 
separation of 50" or ^0.03 pc, which is small compared to 
our estimate for the Jeans length at the two peaks (0.11 
and 0.08 pc for cores 09a and 09b, respectively). Follow- 
ing the prescription of RLA09, the only other possibil- 
ity for the two fragments to represent independent cores 
would be that they are separated in velocity by more 
than the local sound speed. A high resolution C 18 0(2-l) 
map of the central region in B59, recently obtained at the 
IRAM 30m telescope (Roman-Zuniga et al., in prepara- 
tion) allowed us to determine that the central velocities 
at the two peaks barely differ. Our NH 3 single point 
measurements also do not show significant variations of 
the central velocities near the two peaks. Thus, given 
our current definitions, the two peaks at the nucleus of 
B59 would be considered to be part of the same core. 
Moreover, the dust column density distribution is in fact 
consistent with a single, centrally condensed, stratified 
core with two small fluctuations near the top. It is sig- 
nificant that the depression in between the two peaks 
roughly coincides with the location of the protostellar 
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TABLE 2 

CLUMPFIND-2D IDENTIFICATIONS IN HIGH-RESOLUTION EXTINCTION MAP OF BARNARD 59 



Clump ID 1 Peak Ay Sum Ay No. pixels R eq (pc) Mass n Lj Mj 

[mag] in map [pc] [M ] 10 4 cm -3 ] [pc] [M Q ] 



Central Core Region 



09a 


83.3 


11995.8 


1042 


0.11 


10.4 


2.87 


0.12 


1.1 


09b 


83.7 


9796.4 


530 


0.08 


8.5 


6.46 


0.08 


0.7 


09c 


24.0 


3254.0 


554 


0.08 


2.8 


2.00 


0.11 


1.1 


09d 


7.8 


528.7 


209 


0.05 


0.5 


1.41 


0.17 


1.6 


09c 


7.8 


188.6 


105 


0.04 


0.2 


1.41 


0.17 


1.6 



Eastern Filament 



11a 


17.2 


2954.7 


760 


0.10 


2.6 


1.20 


0.18 


1.7 


lib 


14.9 


926.6 


197 


0.05 


0.8 


2.69 


0.12 


1.1 


11c 


14.5 


831.4 


171 


0.05 


0.7 


2.99 


0.12 


1.1 


lid 


11.5 


2017.4 


417 


0.07 


1.7 


1.91 


0.14 


1.4 


11c 


11.0 


1492.9 


347 


0.07 


1.3 


1.86 


0.15 


1.4 


llf 


11.3 


1181.4 


330 


0.07 


1.0 


1.59 


0.16 


1.5 


llg 


7.8 


1741.4 


613 


0.09 


1.5 


0.93 


0.17 


1.9 



Peripheral cores 



03 


19.9 


2747.7 


811 


0.10 


2.4 


0.96 


0.20 


1.9 


04a 


9.4 


887.7 


342 


0.07 


0.8 


1.13 


0.19 


1.8 


04b 


4.6 


778.1 


393 


0.07 


0.7 


0.81 


0.22 


2.1 


05a 


10.9 


804.8 


285 


0.06 


0.7 


1.35 


0.17 


1.6 


05b 


9.9 


1301.0 


396 


0.07 


1.1 


1.32 


0.17 


1.6 


07a 


14.2 


707.5 


259 


0.06 


0.6 


1.36 


0.17 


1.6 


07b 


10.3 


1079.4 


339 


0.07 


0.9 


1.39 


0.17 


1.6 



Core identifications follow the numbering used in the list of RLA09. Identified substructure is 
named after the main RLA09 core number plus letters (a,b,c, etc.). 



Barnard 59. Central Core 



100 



80 



60 - 



20 - 



* 



m 




100 



1000 10000 
Radius from center (AU) ot D p[pe =130 pc 



Fig. 9. — Radial profile of the central core 09ab (data points), 
compared to a King density profile model (thick, dotted curve). 
The center of the core was defined as the mid-point between cores 
09a and 09b. Please see the electronic version of the Astrophysical 
Journal for a full color version of this figure 



source No. 10. This could suggest that the double peak 
near the center of the map could be a consequence of the 
presence of the protostar. 



If 09ab was a single, large core of 19 M© in the pro- 
cess of collapse, it would be expected to produce, at the 
current star formation efficiency of 20%, a star with a 
mass of approximately 4 M©, which would be the most 
massive star produced in the B59 at the current stage of 
evolution. To test the idea that the central core could be 
a single, centrally condensed, quasi-isothermal structure, 
we compared the radial profile of the 09 ab region w ith 
an analytical King density profile model ( |King||l962[ ) of 
the form Y(r) — Yq[(1 + (r /ro) 2 ] -1 , as shown in Figure [9] 
The fit is done in the Ay vs. radius space with units mag 
and AU, respectively. The radial profile of the B59-09ab 
core was constructed with circular, concentric apertures, 
Nyquist sampled and centered in the mid point between 
the two peaks visible at Ay=72 mag. For that reason, 
the profile is depressed near the center. Our model of 
choice fits the core to its peak value, with Yj = 83 mag 
and Tq = 0.85e4 AU. The reduced \ 2 value for the fit 
is 1.81, which represents a moderately good fit, and if 
we compare only the values for r > 1000 AU then the 
reduced \ 2 falls short of 0.97 indicating a good quality 
fit. The agreement of the model with our data is thus 
acceptable, with only the two centermost points deviat- 
ing significantly from the fit, showing a depression. This 
depression near the center of the core, accounts for only 
a 4% difference in mass with respect to the area under 
the model profile. The agreement of the data points with 
the King model reinforces the idea that the core could 
be a single, centrally condensed, nearly isothermal struc- 
ture. At this point, we did not want to compare the 
profile of B59 wit h a Bonnor-Ebert sphere model as in 
Alves et al. (20011. The main reason is that preliminary 
analysis yield a fit with a very large £, m ax , indicative of a 
highly condensed object far from stable equilibrium, and 
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as shown by Kandori et al. (2005), it is complicated to 
distinguish a Bonnor-Ebert sphere in unstable equilib- 
rium from an object in stage of collapse. Such discussion 
has important consequences, but they fall a bit far from 
the scope of this study. In a companion paper (Roman- 
Zuniga et al., in preparation), we compare the profiles of 
cores across the whole Pipe Nebula, including B59, with 
models of singular isothermal spheres, and discuss the 
implications of such analysis. 

5.2. A Quiescent Cluster Forming Core? 

We used our NH3 emission measurements to obtain in- 
formation about the kinematic and thermal properties of 
the dense gas in the center of B59. As expected for a high 
density core, the NH3 emission was found to be relatively 
intense for most of the pointed observations. The GBT 
was pointed towards a total of 37 positions across the 
central core; 26 of these pointed observations yield sig- 
nificant emission in both the (1,1) and (2,2) transitions, 
while the 11 remaining pointings were rejected for having 
a low signal-to-noise ratio in either the (2,2) or the (1,1) 
transition. The 26 final pointings allowed to determine 
velocity dispersions, er„, and kinetic temperatures, T&, 
across the central part of the core via a forward-fitting 
routine (for details on the data reduct i on an d model fit 
methods, please see Rathborne et al. ( 2007 1) . The ki- 
netic temperatures range from 10 to 13 K, with a mean 
of 11.3±0.7, and the velocity dispersions, o~ v varied from 
0.13 to 0.29 with a mean of 0.18±0.05. 

From these measurements we were able to estimate the 
isothermal sound speed, c s and the non-thermal velocity 
dispersion, a n t- Then, we calculated the ratio of thermal 



to non-thermal pressure R„ 



2 

n 1 3 



c s is the isother- 



mal sound speed at the corresponding Tk), and how it 
varies across the center of the core. In Figure [10] we plot 
the ratio R p vs the extinction value at the center of the 
pointing. Errors were estimated as follows: in the case 
of R p we took into consideration the uncertainties of the 
forward fitting routine for Tk and a v ; in the case of Ay 
we considered both the value of the median noise at the 
value of the pixel closest to each and the mean absolute 
deviation of the Ay values in the 9 surrounding pixels 
to account for the uncertainty due to the differences in 
beam sizes (GBT:30" vs. NICER map:24"). 

It is certainly striking that a majority of the values 
are greater than R p =l, indicating that thermal motions 
dominate over non-thermal (turbulent) motions of the 
gas in the central core. The mean of R p in B59 is 
1.90±0.99; for comparison, a thermally dominated core 
like FeSt-1457 (which also belongs to the Pipe Nebula), 
was shown to hav e an average ratio of R p = 3.75 ± 1.95 
( |Aguti et al.|2007[ ). The plot of Figure [lO] also shows that 



-with ameanof R p — 1.0±0.4- 



the 12 lowest values of R. 

are at the densest part of the core (log Ay >1.5 mag). 
This suggests that non-thermal motions make a more sig- 
nificant contribution to the internal pressure at the cen- 
tral part of the core, where the column density is higher. 
The relatively large R p ratios in B59 tell us that the cen- 
tral core is still thermally dominated despite the pres- 
ence of 14 cluster members packed in the same observed 
region. Moreover, most of the non-thermal velocity dis- 
persion values in B59 are smaller than the correspond- 
ing isothermal sound speed, suggesting that this cluster 
forming core is characterized by sub-sonic turbulent mo- 
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Fig. 10. — Top panel: average spectra of NH3 (1,1) and (2,2) 
emission towards the center of B59. Bottom panel: Thermal to 
non-thermal pressure ratio vs. visual extinction in B59. 



tions. 

5.3. Implications 

Our observations appear to indicate that a) there is 
a lack of significant sub-structure in the central region 
of B59 and b) that the central core is mostly quiescent, 
and still thermally dominated. The result is somewhat 
unanticipated, as we expected to find evidence of mul- 
tiple fragmentation in a core with M ~ 20Mj. More- 
over, the significant sub-structure we do detect in the 
central core of the B59 system appears to be related to 
the embedded protostars, while the overall radial pro- 
file of the core seems to be consistent with one large, 
isothermal structure. Therefore, the question that arises 
is: what prevents B59 from fragmenting and collapsing 
into new stars after a period considerably larger than its 
own dynamical timescale? In one theoretical scenario, 
a first generation of stars is formed, and they provide 
feedback which would contribute to suppress fragmenta- 
tion in the remaining gas, allowing for accret ion of ma- 
terial to form a more massive central object fMcKee fc| 
Tan||2002| [krumholz k, I an||2007|. Recent observations 
show that the addition of significant non-thermal mo- 
tions could prevent fragmentation and allow for the ag- 
gregation of larger amounts of material that co uld later 
collapse into a massive star ( |Zhang et al.||2009| . 

As we discussed before, our map shows three possi- 
ble signatures of stellar feedback in B59: the first is the 
large northwestern cavity being created by source No. 11, 
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which might have triggered the formation of core 09c; the 
second is the other cavity-like indentation at the south- 
east edge of the central core, possibly associated with 
source No. 9; the third is the depression at the center 
of the core, associated with source No. 10. However, we 
cannot detect significant non-thermal motions produced 
by the stellar feedback in the linewidths of the radio ob- 
servations. Therefore, additional non-thermal support 
from some other process might be required. 

One possibility is that additional support against col- 
lapse could be provided by a mech anism of ambipolar dif- 
fusion via a strong magnetic field. Kauffmann & Bcrtoldi 
( 2004 1 showed that massive molecular cores tend to have 



larger column densities than expected from pure thermal 
support, but also that in some cases the non-thermal 
linewidths are smaller than expected (as in B59), in 
which case, support from a strong magnetic field would 
be required. The hypothesis of a strong magnetic field 
acting in the Pipe Nebula has been supp orted by nu- 
merical simulations ( Nakamura fc Li||2008 l, which show 
that ambipolar diffusion is an agent capable of slowing 
down the process of condensation, as needed in the Pipe, 
where the free-fall time scale is very short given the me- 
dian core density. A magnetic field could provide the 
necessary non-thermal support against collapse in mas- 
sive cores like B59, giving enough time for turbulence to 
decay and allowing the existence of quiescent condensa- 
tions with moderate supersonic motions, similar to what 
we observe. A recent R-band o ptical polarimetry study 
of the Pipe by Alves et al. ( 2008 1 suggested that the mass 
to magnetic flux ratio toward 1359 is about 1.4 (slightly 
supercritical), with a field strength 2 to 3 times weaker 
than other parts of the cloud. However, note that most of 
the stars available for optical polarimetry in B59 are lo- 
cated in regions of low column density at the edge of the 
core, and thus a more reliable tracer of the field towards 
the core could still improve these estimations. 

5.4. Fragmentation the Eastern Filament 

In the RLA09 analysis of the 2MASS extinction map, 
the eastern filament in B59 was considered to be one 
large core because their measurement of the length and 
separation of the fragments detected by CLUMPFIND-2D 
did not comply with the requirements of being separated 
by more than a Jeans length and by more than 0.12 km/s 
in central velocity. Our high resolution map allows us to 
revisit the first requirement with a better quality dataset. 

In Figure [IT] we show the positions of sources in the 
area of the filament with reliable quality for application 
of the near-IR color excess method, from the 2MASS 
catalog and the ESO-NTT survey observations. The dif- 
ference in the density of sources is remarkable, but not 
as much as the fact that the central region along the fil- 
ament is almost void of sources in 2MASS. This means 
that LAL06 were probably not able to resolve the cen- 
tral part of the filament. Instead, with the ESO-NTT 
data we are sampling extinction all the way to the mid- 
dle region. While the 2MASS map and the ESO-NTT 
map agree in the basic morphology, extension and basic 
structure of the core, the 2MASS map can not be used 
to det ermine the presence of substructure: |Lada et al.| 
( |1999| ) showed (specifically for a cylindrical cloud), that 
a poor sampling of the background population in a highly 
obscured region due to a limited photometric, depth can 
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Fig. 11. — A comparison of the density of background sources 
from the 2MASS (top) and the ESO (bottom) catalogs, in the 
region of the eastern filament. 



lead to the creation of spurious sub-structure, or at least 
to an unreliable determination of the size and mass of 
substructures. In that sense, our hybrid map resolves 
the filament along its axis and allows for the recognition 
of significant substructure. 

Our analysis of the filament with CLUMPFIND-2D re- 
veals a total of 6 extinction peaks along the filament, 
and that these peaks are significant with our choice of 
5 sigma contour steps. However, we must proceed with 
caution: again, it not trivial to claim that each one of 
these peaks should be considered a new core. The RLA09 
prescription requires individual cores to be separated in 
both distance and velocity. In their analysis, RLA09 
also detected substructure along the filament, but they 
decided to merge all peaks into a single, large core as 
the adjacent peaks in the map of LAL06 did not show a 
separation in central velocity However, we need to con- 
sider that while the prescription of Rathborne is useful 
to determine the independence of two adjacent cores, it 
might not be ideal for a chain of related cores like in this 
filament. 

We suggest that the filament might be in a process 
of thermal fragmentation: we do not have pointed radio 
measurements at the peaks in the high resolution map, 
but we can revise the condition of the projected distance 
between the peaks along the filament with our improved 
measurements of the LJL: the small projected separation 
between peaks lib and 11c is smaller than the LJL, sug- 
gesting that they might represent a single fragment, and 
thus we merged them in the final list of cores. Our peak 
to peak measured distances for the remaining features 
average 0.11±0.3 pc, just shorter the local Jeans length 
estimated for the individual peaks, which varies from 0.12 
to 0.16 pc. We have to consider, however, that projection 
effects over such a long structure may not be small and 
thus our measured projected distances are lower limits. 
Also, the resultant Jeans mass values, Mj for the peaks 
along the filament average to 1.3±0.2 M , and the ra- 
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tio Mj/M has an average of 1.2±0.1, indicating that the 
peaks along the filament are Jeans stable. Moreover, we 
measured the projected width of the filament directly on 
the map at seven different points along the eastern fila- 
ment. The average width is 0.10±0.02 pc, that compares 
well with the average LJL across the filament. 

These measurements suggest that the eastern filament 
in B59 could be following a process of thermal fragmen- 
tation, and could break into several fragments. The 
fragmentation of isothermal, cylindrical clouds has been 
studied often: observations by |Schneider fc Elmegreen 



(19791 suggested that in a majority of observable fila- 



ments in molecular clouds, it was possible to observe a 
periodicity in the separations of peaks along filaments, 
and that in many cases these fragments agreed spa- 
tially with the locations of young sources, suggesting that 
peaks observed along filaments are the precursors of in- 
dependent cores. Recent observations have shown signif- 
icant eviden ce of thermal fragment a tion of filamentary 
clouds (e.g. |Wiseman fcH oj jl998j); |Hartmann| ( |2002| ); 
Teixeira et al.| ( |2005[ |20rj5| )~^ 

Filamentary structures may become unstable and col- 
lapse via thermal fragmentat ion once they lose equilib- 
rium with external pressure. McCrea (1957) found that 



Bastienj(|1983| allow 
the Jeans number, 



a finite cylinder would collapse once its mass per unit 
length, m = M/L surpasses a critical value, m c . Numer- 
ical two-dimensional simulations by 
one to make a rough estimation or 
J cy i, in a finite isothermal cylinder of length L and di- 
ameter D as J cy i « 3(m/m c )f(L/D), where f(L/D) is 
a geometrical factor. Following Bastien calculations, for 
a ratio L/D w 5 as observed in our map, the value of 
f(L/D) from Bastien tables is close to 1.6, the m/m c ra- 
tio is about 0.3, and the critical Jeans number is J c =0.6, 
so that J cy i I J c = 2.4. The expected number of fragments 
is then roughly equal to (J cy i/J c ) 3 ^ 2 ~ 3.6, not far from 
our estimate of five fragments along the eastern filament. 

5.5. A revised list of cores in Barnard 59 

In table [3] we re-organize the list of map features ac- 
cording to the previous discussion. The table also list the 
locations of the cores and their peak extinction values 
in the unfiltered map. We confirm the seven main core 
identifications made by RLA09 in the map of LAL06 The 
analysis of our high resolution map confirms that cores 
03 and llg are single globules, while cores 04, 05, 07, 09 
and 11 have significant substructure. In the case of core 
11, the eastern filament, we find five fragments of sig- 
nificance. We suggest that peaks lib and 11c represent 
one single fragment and that core llg might not belong 
to the same chain. Regarding the central core we con- 
clude that it has little observable fragmentation despite 
its large mass. The central core has two significant frag- 
ments, 09ab and 09c, the latter possibly being formed by 
a protostellar outflow; cores 09d and e are merged into a 
single independent core, 09de. 

6. SUMMARY 

The main results we achieved in this study can be listed 
as follows: 

1. Using a combination of near-infrared and mid- 
infrared (3 to 8 yitm) observations, we obtained very deep 
photometry for the dense population of background stars 
(mostly galactic bulge giants) towards Barnard 59, the 



most dense and dark region of the Pipe Nebula. The 
dense central core of the Barnard 59 complex is one of 
the very few cores in the cloud with active star formation 
and the only one forming a stellar cluster. 

2. Our deep observations allowed us to make a de- 
tailed map of dust extinction towards the B59 region 
with unprecedented spatial resolution (12 to 24") and 
dynamic range of column density (Ay < 89 mag or 
N H2 < 8.45 x 10 22 cm- 2 ). 

3. Our extinction map allowed us to determine, in de- 
tail, the structure of the central core as well as of the 
filaments and peripheral cores in the Barnard 59 com- 
plex. We found a total of 13 independent extinction 
peaks in the region, most of them interconnected by a 
complex structure of filaments. The largest filament is 
characterized by a high column density, with a peak ex- 
tinction above Ay=22 mag. It appears to contain five 
independent cores evenly spaced along its axis. 

4. The total mass of the complex is approximately 41 
Mq, with half of this mass (about 21 M Q ) being con- 
tained in the central region (comprising cores 09ab and 
09c), which is associated with the formation of the stellar 
cluster. 

5. The central core appears to have remarkably little 
structure despite the fact that its mass is significantly 
larger than either the Bonnor-Ebert or the Jeans mass 
in the region. The split of the central core into two peaks 
is significant, but the small projected separation of the 
peaks and the negligible difference in central velocities 
does not support the idea that these peaks represent two 
independent cores. Alternatively, the peaks could be re- 
lated to the presence of a class /I source which is located 
between them. This source may have created a depres- 
sion in what otherwise would be a large stratified core 
with a 1/r 2 density distribution. The apparent separa- 
tion of a fragment at the northwestern edge of the central 
core also appears to be significant, but it is at the edge of 
an outflow cavity and could be the result of the disrup- 
tion of the core material by the stellar outflow. Our map 
also suggests that a second cavity is located near source 
No. 9. Apparently, most of the observed substructure 
in B59 could be the result of action by embedded YSOs, 
indicating that stellar feedback is playing role of signifi- 
cance in the evolution of the core. 

7. Pointed observations of the (1,1) and (2,2) rota- 
tional transitions of NH3 at 26 positions in the central 
core show that it is thermally dominated, with non- 
thermal motions being mostly subsonic. The central 
core has formed a small cluster of stars, but presently 
it is mostly quiescent and does not show any kinematic 
evidence of global collapse. It is possible that the feed- 
back from the recently formed stars may be preventing 
the global collapse and suppressing further fragmenta- 
tion. However, the absence of a significant non-thermal 
component in the observed radio lines suggests that an 
additional source of non-thermal support is required. 

8. Substructure along the axis of the eastern filament 
appears to be consistent with Jeans fragmentation. 
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TABLE 3 

Revised List of Cores in Barnard 59 



Clump ID 1 RA Dec Peak ExtinctionID 2 C ore Mass 

[J 2000] [mag] [M Q ] 



Central Core Region 



09ab 


17:11:21.8 


-27:26:10.0 


89.0 


18.9 


09c 


17:11:27.1 


-27:23:10.0 


29.8 


2.8 






Eastern Filament 






11a 


17:12:32.2 


-27:21:31.8 


20.0 


2.6 


llbc 


17:12:21.9 


-27:23:49.6 


18.3 


1.5 


lid 


17:12:09.3 


-27:25:19.8 


15.7 


1.7 


lie 


17:11:58.8 


-27:26:22.4 


16.2 


1.3 


llf 


17:11:50.3 


-27:25:02.5 


16.2 


1.0 






Peripheral Cores 






03 


17:10:31.5 


-27:25:46.6 


22.4 


2.4 


04 


17:11:35.9 


-27:33:30.0 


13.4 


1.5 


05 


17:12:13.7 


-27:37:19.7 


13.1 


1.8 


07 


17:10:51.8 


-27:22:59.6 


13.3 


1.5 


09de 


17:11:14.1 


-27:22:32.4 


12.8 


0.7 


Hg 


17:12:52.1 


-27:23:08.7 


11.0 


1.5 



Core ID after merging of features as described in text 
2 Refers to peak extinction of core in the unfiltered map 
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